Graphene on metals: a Van der Waals density functional study 
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We use density functional theory (DFT) with a recently developed van der Waals density func- 
tional (vdW-DF) to study the adsorption of graphene on Al, Cu, Ag, Au, Pt, Pd, Co and Ni(lll) 
surfaces. In constrast to the local density approximation (LDA) which predicts relatively strong 
binding for Ni,Co and Pd, the vdW-DF predicts weak binding for all metals and metal-graphene 
distances in the range 3.40-3.72 A. At these distances the graphene bandstructure as calculated with 
DFT and the many-body Go Wo method is basically unaffected by the substrate, in particular there 
is no opening of a band gap at the if-point. 

PACS numbers: 71.15.Mb,71.15.Nc,73.20.Hb 



The recently reported synthesis of graphene [lj], a sin- 
gle layer of graphite, ontop of a Si02 substrate has re- 
newed the interest for this unique material. The unique- 
ness of this 2D crystal is mainly due to its very peculiar 
band structure, with the n and ir* bands showing linear 
dispersion around the Fermi level where they touch in 
a single point. The great variety of physics and chem- 
istry which derives from this electronic structure makes 
graphene very attractive for a range of applications. In 
particular, its high stability and good conductivity un- 
der ambient conditions makes it an interesting candi- 
date for future nano-scale electronics In this perspec- 
tive, the interaction of graphene with metallic contacts 
plays a fundamental role. Moreover, catalytic growth of 
graphene on transition metal surfaces from carbon con- 
taining gases has become a standard way to obtain high 
quality graphene samples Nevertheless the na- 

ture of the metal-graphene chemical bond is still not well 
understood. 3 . 

The widely used density functional theory (DFT) with 
local and semi-local functionals for exchange and correla- 
tion usually provides an accurate description of covalent 
and ionic chemical bonds. On the other hand it fails to 
reproduce non-local dispersive forces, in particular van 
der Waals forces, which are important in weakly bonded 
materials such as graphite, molecular crystals, and many 
organic compounds. 0, [1, @|. It is also well known that 
the local density approximation (LDA) tends to overbind 
systems where van der Waals interactions are important, 
while the generalized gradient approximations (GGA) 
usually tend to underestimate the binding in these sys- 
tems. In the case of graphene on metals many GGAs, 
contrary to experiments, predicts no binding at all, and 
therefore most theoretical work on graphene-metal inter- 
faces has relied on the LDA. In view of the fact that 
LDA in general cannot be considered a reliable approxi- 
mation in non-homogeneous systems such as surfaces and 
molecules, the graphene-metal interface clearly calls for 
new and improved functionals. 

The interaction of graphene with the (111) surfaces of 



Al, Cu, Ag, Au, Pt, Pd, Co and Ni was studied in Ref.[10( 
using the LDA approximation. The LDA results divide 
the metals into two classes: Co, Ni and Pd which bind 
graphene strongly and Al, Cu, Ag, Au and Pt which bind 
graphene weakly. In contrast PB Ella gives no binding 
of graphene at room temperature [l3(. This remarkable 
disagreement between the two most commonly used ap- 
proximations of DFT might be related to the incorrect 
description of dispersion interactions in both of the func- 
tionals. 

In this paper we use the recently dev elop ed van der 
Waals density functional (vdw-DF)[3, EEf to investi- 
gate the nature of the bonding at the metal-graphene 
interface. The functional is explicitly constructed to in- 
clude non-local dispersion interactions and has proven 
successful in several cases where standard functionals fail, 
such as rare gases 14 [, benzene dimersfltj 17 1, graphite 
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the exchange-correlation energy is 
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where ££= vPBE is the revPBEfll] exchange energy, £^ DA 
is the LDA correlation energy and E£ is the non-local 
correction given by 



Ef = i J J n(r 1 )n(r 2 )(f>(q 1 ,q 2 ,ri2)dr 1 <ir 2 



(2) 



where r\i = \x\ — Vi\ and q\ and q 2 are values of a uni- 
versal function q (n(r), | Vn(r)|). Eq. ^ is efficiently 
evaluated by factorizing the integration kernel <p and by 
using fast Fourier transform to compute the selfconsistent 
potential as proposed in Ref. [13] and implemented in the 
real-space projector augmented wave GPAW code[30]. 

In this study we consider Al, Cu, Ag, Au, Pt, Ni, Co 
and Pd metal (111) surfaces. We fix the atoms in the 
metal slabs at their experimental lattice parameters and 
relax the graphene sheet using the vdw-DF Hellmann- 
Feynman forces. We use a (6,6,1) and (4,4,1) Monkhorst 
Pack k-point sampling respectively for the smaller (Ni, 
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TABLE I: Binding energies (Eb) per carbon atom and binding 
distances (d) of graphene on metal (111) surfaces; '-' means 
no binding. Fermi level shift AEf and charge transfer 8Q 
of graphene adsorbed on the different metals at the vdw-DF 
equilibrium separation. Negative (positive) AEp indicates n 
(p)-type doping. Negative (positive) SQ indicates electron 
transfer to (from) the graphene layer. The char ge t ransfer 
has been evaluated according to the Bader scheme [311] . 



Cu, Co) and the larger (Au, Ag, Pt, Pd, Al) orthorombic 
unit cells. The metal slabs are modeled with 4 atomic 
layers and a vacuum of 14 A in the direction normal to 
the surface; the grid spacing is 0.16 A. The calculated 
binding energies and distances for the relaxed structures 
are listed in Table [TJ The vdW-DF results show that the 
metal-graphene interaction is relatively similar across the 
different metals. This is in contrast with the LDA pre- 
diction of two separate classes of metal-graphene inter- 
faces, as found in very good agreement with (Toj . We 
also repeated the same calculations using the revPBE 
functional and we obtained no binding for any of the 
metals. Interestingly, for the systems that LDA finds to 
be weakly bonded (Al, Cu, Ag, Au and Pt), the binding 
energies obtained with the vdw-DF are very similar to 
the LDA ones. Nevertheless the binding distances are 
systematically slightly larger in the vdw-DF case. In 
fact it has been reported that the vdw-DF functional 
usually produces equilibrium distances somewhat larger 
than experiments 2l|. In the case of Ni, Co and Pd, on 
the other hand, the relatively strong binding predicted 
by LDA is not found by the vdw-DF functional. 

In order to analyze these results, we now focus on the 
interaction between graphene and Ni(lll). Fig. [1] shows 
the binding curves for graphene on the Ni(lll) surface 
calculated with the LDA, revPBE and vdw-DF function- 
als. The revPBE curve is positive at all distances, while 
the LDA curve shows a relatively deep minimum at ~ 2 A 
consistent with previous LDA calculations. The vdw- 
DF result lies in between, following the revPBE curve at 
small separations and the LDA curve at larger separa- 
tions, and it predicts a shallow minimum at 3.5 A. Note 
that a local minimum is found by the revPBE functional 
around 2 A. 

In Fig. [2] we show the calculated bandstructure of 
graphene on Ni(lll). The color of the dots indicate the 
weight of the corresponding Bloch eigenstate on the car- 
bon p z orbitals with darker meaning larger weight. In 




FIG. 1: Binding energy (Eb) per carbon atom of graphene on 
the Ni(lll) surface calculated with LDA, revPBE and vdw- 
DF functionals. The graphene is adsorbed in the top-fcc con- 
figuration. 



free graphene, the carbon p z orbitals placed at A sites 
(p A ) are decoupled from the p z orbitals at B sites (pf ) 
at the Dirac point, thus producing two degenerate states 
(see inset in Fig. [T]for the structure). Since the A sites 
are located directly on top of Ni atoms at a close distance 
in the LDA calculation (2.08 A), a strong hybridization 
between p^ orbitals and Ni 3z 2_ r 2 is observed, which gives 
rise to an unoccupied antibonding state a* and two occu- 
pied bonding states a± and . The LDA gaps for a* — o\ 
and a* — o"2 are 2 eV and 4 eV, respectively. On the other 
hand, the p^ orbitals (occupied in the spin up channel 
and unoccupied in the spin down one) hardly interact 
with Ni d states and therefore remain unmodified. The 
vdw-DF band structures (evaluated at the vdw-DF re- 
laxed distance of 3.50 A), on the other hand, resemble 
the free graphene, preserving the Dirac point and only 
shifting it up by 0.13 eV. A very similar behaviour is 
found for Co and Pd. For the remaining interfaces both 
the LDA and vdw-DF bandstructures resemble that of 
free standing graphene with the Dirac point shifted with 
respect to the metal Fermi level. The Fermi level shifts 
and calculated charge transfer between the metal and 
graphene are summarized in Table [I] 

Since LDA is known to underestimate band gaps signif- 
icantly we have also performed Go Wo calculations for the 
graphene-Ni structures corresponding to the LDA and 
vdw-DF distances. [25j In both cases we find no notewor- 
thy difference between the Go Wo and DFT Kohn-Sham 
band structures close to the ii-point. 

Recent experimental work on the graphene/Ni inter- 
face is based on Angle Resolved Photoemission Spec- 
troscopy (ARPES). The ARPES band structures re- 
veal a band gap at the graphene ii-point thus suggest- 
ing some h ybr idization between the graphene and Ni 
orbitals 27, 281. Earlier LEED measurements found a 
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FIG. 2: LDA (top) and vdw-DF (bottom) spin polarized band structures for graphene on Ni(lll) in the top-fee configuration. 
Darker dots represent larger weight of the carbon p z orbitals. 




FIG. 3: Decomposition of the binding energy Ef, into exchange-only contributions (dashed lines - only the correlation term is 
removed) and total binding energy (full lines) for different functionals. 



Ni-graphene bond distance of 2.lA[29(. We note that 
both of these results are in line with the LDA calcula- 
tions. On the other hand LDA is not expected to work 
well for highly inhomogeneous systems such as the inter- 
face structures investigated here. These results indicate 



a difficulty for the present vdw-DF in describing systems 
with mixed bonding character, in line with the conclu- 
sions of [23j |. 

Fig. [3] shows the total (full lines) and the exchange- 
only (dashed lines) binding energy curves for revPBE, 
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PBE and LDA in the case of graphene on Cu(lll) (left 
panel) and on Ni(lll) (right panel). The exchange-only 
energies are calculated without including the correlation 
energy term in the exchange-correlation functional and 
have been evaluated non selfconsistently. Clearly, the 
bonding for the physisorbed graphene on Al, Cu, Ag, Au 
and Pt originates partially from the exchange term in 
the LDA xc-functional, as shown for Cu as an example 
in the left panel of Fig. [3] . This is in principle incor- 
rect since van der Waals interactions are purely non-local 
correlation effects. The weak bonding predicted by the 
vdw-DF functional, similar in magnitude to the LDA re- 
sults, is produced by the correlation term instead, which 
is physically correct. Interestingly, this applies to the 
Ni/graphene system as well, as shown in the right panel 
of Fig. [3J. The LDA exchange-only curve shows a broad 
and weak attractive contribution between 2 A and 5 A 
with two local minima. We note that in a genuine co- 
valent bond the exchange contribution to the binding 
energy is generally significantly larger than in this case. 
The revPBE and PBE exchange-only curves are repulsive 
at all separations for both systems. This is the typical 
behaviour which is observed in van der Waals bonded 
dimers or organic molecules on surfaces 3, 22 1 

In conclusion we have performed DFT calculations of 
graphene adsorbed on different metal surfaces using the 
recently developed vdw-DF functional which explicitly 
includes non-local correlations. For Al, Cu, Ag, Au and 
Pt both LDA and vdw-DF consistently predicts a weak 
binding. Similar weak binding is found for Ni, Co and 
Pd with the vdw-DF where LDA on the other hand 
predicts stronger binding and significant hybridization 
between graphene and metal d-states. At the vdw-DF 
binding distances graphene's band structure was shown 
to be essentially unaffected by the substrate. This ap- 
pears to be in conflict with LEED and ARPES measure- 
ments for graphene on Ni, indicating that more work is 
needed in order to reconcile experiments and theory for 
the graphene-metal interface. 
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